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ABSTRACT 



Context. We investigate the deposition of energy and momentum due to the annihilation of neutrinos (i^) and antineu- 
trinos in the vicinity of steady, axisymmetric accretion tori around stellar-mass black holes (BHs). This process is 
widely considered as an energy source for driving ultrarelativistic outflows with the potential to produce gamma-ray 
bursts. 

Aims. We analyze the influence of general relativistic (GR) effects in combination with different neutrinosphere prop- 
erties on the i/j5-annihilation efficiency and spatial distribution of the energy deposition rate. 

Methods. Assuming axial symmetry, we numerically compute the annihilation rate 4-vector. For this purpose we con- 
struct the local neutrino distribution by ray-tracing neutrino trajectories in a Kerr space-time using null geodesies. We 
vary the value of the dimensionless specific angular momentum a of the central BH, which provides the gravitational 
field in our models. We also study difi'erent shapes of the neutrinospheres, spheres, thin disks, and thick accretion tori, 
whose structure ranges from idealized tori to equilibrium non-selfgravitating matter distributions. Furthermore, we 
compute Newtonian models where the influence of the gravitational field on the annihilation process is neglected. 
Results. Compared to Newtonian calculations, GR effects increase the total annihilation rate measured by an observer 
at infinity by a factor of two when the neutrinosphere is a thin disk, but the increase is only « 25% for toroidal and 
spherical neutrinospheres. Comparing cases with similar luminosities, thin disk models yield the highest energy depo- 
sition rates by i^i/-annihilation, and spherical neutrinospheres the lowest ones, independent of whether GR effects are 
included. Increasing a from to 1 enhances the energy deposition rate measured by an observer at infinity by roughly a 
factor of 2 due to the change of the inner radius of the neutrinosphere. General relativity and rotation cause important 
differences in the spatial distribution of the energy deposition rate by !/i?-annihilation. 

Key words, gamma-rays; bursts - neutrinos - accretion, accretion disks - relativity - black hole physics - stars: neutron 



1. Introduction 

It is widely believed that systems powering gamma-ray 
bursts (GRB) could be newborn, stellar-mass black holes 
(BHs) accreting matter at hyper-critical rates (up to sev- 
eral solar masses per second) from a surrounding accretion 
disk with a mass of some hund redth of a so lar mass up to 
possibly a solar mass (see, e.g.. lPiranll2005[ ). These central 
engines may form in a 'coUapsar' event where the core of a 
massive, rotating Wolf-Rayet star collapses to a BH and the 
accretion of the stellar envelope may eventually lead to a 
GRB-supernova ev ent with relativ i stic mass ejection alon; 
the rotation axis (IWqoslevI [19931 iMacFadven fc Woosle 
Il999t lAlov et all 120001). Accreting BHs may also be the 
remn ants of mergers of two compact objects in clo se bina- 
ries (jEichler et al.lll989t [Mochkovitch et~allll993[ ). In the 
first scenario the system is embedded in the envelope of 
the progenitor star, and the accretion disk is fed by stellar 
matter yielding very long (^ 10 — 1000 s) accretion time 
scales comparable to the collapse time scale of the progen- 
itor star. In the second scenario viscous transport in the 
accretion torus sets the secular time scale of the system 
(~ 0.01 — Is) formed during the merger. 

The conditions in the vicinity of steady-sta te, hyperac- 
creting BHs have been analytically studied by I Jaroszvnskil 



(119931. ll996f);lPqp ham et all (Il999l): iNa ravan et al.' (2mh; 
Pi Matteo et all ( 2002): iKohri fcl-Iincshigc (200^; 
Chen fc Beloborodovi p006[ ). who determined the ef- 
ficiency of energy loss by neutrino emission and the 
efficiency of energy conversion by neutrino-antineutrino 
(I'l') pair annihilation into electrons and positrons. Three- 
dimensional hydrodynamic simulations have explored the 
time- dependent accretion in BH-torus systems, which are 
the remnants of n eutron star-neutron sta r (NS-I-NS) and 
NS-f- B H mergers (iRuffert fc Jankal 119991: ISetiawan et"all 



I2004t iLee et al.l l2004 l2005al lbr). These 



investigations 

considered a rather compact torus (typical size: ^10—20 
Schwarzschild radii) containing between a few hundredth 
of a solar mass and some 0.1 M©. Such torus masses 
result from NS-I-NS and NS-I-BH merger simulations 
(iJanka et al.lll999l: iRuffert fc J anka"l9 99l: iJanka fc RuffertI 



20021: iRosswog et all 



20031: IShibata et all l2003l . 12005 



Oechslin fc Jankal 120061 )" The torus is partly opaque to 



neutrinos because of its high density. Typically, neutrino 

Under 

ly ^ C + c" ^7 + 7 
in the close vicinity of 



luminosities in excess of 10'^" erg s 
these conditions, the reactions v 4 
give rise to energy deposition 



are produced. 
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the BH at rates rangi ng from several 10 ergs up to 
more th an 10^^ ergs" ^ ("Ruffcr ^ Jankall 19991: IJanka et all 
[1999; S etiawan et al.i r2004. 200^. The resulting e+e'-pair 
plasma-photon fireball may power an ultrarelativistic out- 
flow of baryons with typical Lorentz factors of 10^ — 10'^, 
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provided the baryon loading (i.e., the baryon rest mass 
compared to the internal energy in t he outflow) remains 
sufBcicntly low (see, e.g., lAlov et all [2b05i and references 
therein). 

Several studies involving different levels of sophistica- 
tion have been performed, in order to determine the amount 
of energy which can be released by i/z/-annihilation near a 
(ro tating) stella r -mass BH. 

iJaroszvnskil (|l993l I1996D investigated stationary con- 
figurations consisting of massive, dense and non- 
selfgravitating tori with different angular momentum dis- 
tributions and different specific entropies, orbiting stellar- 
mass Kerr BHs. In particular, he considered the neutrino 
emission of isentropic tori in the external field of Kerr 
BHs as an approximation of what results from merger 
ev ents or failed supernovae. Using the neutrino opacities 
of iBurrows fc Lattimeil (11986*) and assuming that the neu- 
trinos have an equilibrium (Fermi-Dirac) distribution given 
by the te mperature and che mical potential at the neutri- 
nosphere, IJaroszvnskil (|1996f ) determined the neutrino ra- 
diation field at a given point by following backwards null 
geodesies in Kerr spacetime until a point is reached in the 
torus where the neutrino optical depth is r^/ « 1. He found 
that the energy deposition rate due to i/C'-annihilation in- 
creases with increasing entropy and with increasing spin 
of the Kerr BH. He also claimed that BH-torus configura- 
tions resulting from NS-I-NS merger events do not provide 
sufficient energy to likely be sources of GRBs. 

Extending e ar lier Newto nian calculati ons of 

Cooperstein et all (I1986D and iGoodman et all (|1987D . 



Salmonson fc WilsonI (|l999t ) analytically determined the 



proper energy deposition rate per unit proper time by 
i^P-annihilation near the surface of a spherical NS (the 
neutrinosphere was assumed to coincide with the NS 
surface) including general relativistic (GR) effects. They 
concluded that the inclusion of the GR effects of ray bend- 
ing and redshift (neutrino trajectories were computed in 
the Schwarzschild metric) enhances the proper deposition 
rate per unit of proper time compared to the Newtonian 
values by up to a factor of 4 for a neutrinosphere located 
at 2.5 Schwarzschild radii relevant for a proto-NS. An 
enhancement factor of 30 is possible, if the radius of the 
neutrinosphere shrinks to 1.5 Schwarzschild radii, which 
ha ppens during the collapse o f a NS to a BH. 

lAsano fc Fukuvainal (|2000t ) studied the influence of GR 
effects on the i/i'-annihilation rate assuming two differ- 
ent geometries of the i^iz-emitting regions: a spherically 
symmetric emission region, and a thin disk emitter sur- 
rounding a Schwarzsch ild BH. The spheric a l geo metry 
was already studied bv I Salmonson fc WilsonI ()1999f ). but 
lAsano fc Fukuvamal (l200(F improved on this work by ap- 
proximately taking into account that some fraction of the 
deposited energy might not escape from the gravitational 
potential well, and thus cannot contribute to power GRBs. 
For disk-shaped neutrinospheres they determined the en- 
ergy deposition rate per unit world time for an observer 
located at infinity by computing the annihilation rate near 
the symmetry axis. The latter restriction allows for a semi- 
analytic treatment by making use o f the axial symmetry 
of the neutrino source. Contrary to ISalmonson fc WilsonI 
()1999f ). they found that GR effects do not substantially 
change the energy deposition rate, neither for the spher- 
ically symmetric case nor for the disk case. This discrep- 
ancy arises, because the energy deposition rates were cal- 



culated with differen t cons traints in both investigations. 
ISalmonson fc WilsonI ()1999D considered the proper energy 
deposition rate per unit proper time, which is enhanced 
by the GR effects of gravitational blueshift and bending 
of trajec tories, for a fixed n e utrino luminosity at infin- 
ity, while lAsano fc Fukuvamal (|2000D computed the (local) 
energy deposition rate per unit world time by assuming 
a given value of the effective temperature of the neutri- 
nosphere and thus of the neutrinosphere luminosity. The 
inverse redshift factor relating the local luminosity to the 
luminosit y at infinity explains the enh ancement of the rate 
m case oflSalmonson fc WilsonI (Il999[). However, as argued 
bv lAsano fc Fukuvama (2000 ). the neutrinosphere luminos- 
ity, which is restricted or provided by models of the GRB 
engine, is the appropriate quantity to use when determining 
the influence of GR effects on the z/i?-annihilation process. 

Im proving on their earlier work, lAsano fc Fukuvamal 
(|200lD considered a non-selfgravitating, geometrically in- 
finitely thin accretion disk surrounding a Kerr BH. The 
disk was either assumed to be isothermal or to have a pre- 
scribed temperature gradient. They found that if the en- 
ergy deposition (along the rotation axis) is mainly due to 
neutrinos from the central part of the disk near the hori- 
zon, redshift effects are dominant, and the energy depo- 
sition rate is consequently reduced compared to the case, 
in which GR effects were neglected during the computa- 
tion. Instead, if neutrinos that are emitted at larger radii 
dominate the energy deposition rate, GR bending effects 
become important. This causes a GR enhancement of the 
energy deposition rate by a factor of about 2, irrespective 
of the spin of the Kerr BH. 

iMiller et all ^qM) considered the full three- 
dimensional problem of numerically computing the 
I'P-annihilation around a Kerr BH / thin accretion 
disk system using the full geodesic equations. They 
determined the energy-momentum deposition rate 4- 
vector per unit 4-volume in a given local observer's 
orthonormal frame. From this covariant quantity they 
derived the (coordinate-dependent, i.e. nonconvariant) 
energy-momentum deposition rate per proper time in 
Boyer-Lindquist coordinates for their observer, which is an 
adequate approximation to the 4-momentum per proper 
time at the observer's location as lon g as the observer is 
not too close to the horizon of the BH jMiUer etall l|2003[) 
imaged the accretion disk for a specified, but arbitrary 
off-axis observer. They did this by integrating geodesies 
along the Boyer-Lindquist energy-momentum vectors until 
they either hit the BH, the disk, or reach r = 5QM. 
They confirmed the findings of lAsano fc Fukuvamal (|2001f ) 
as to an only moderate GR enhancement of the energy 
deposition rate near the rotation axis, but they also showed 
that the dominant contribution to the energy deposition 
rate comes from near the surface of the disk, where the 
rate is a factor 10-20 times larger than on the rotation 
axis. This enhance ment of the r ate i s independent of 
the spin of the BH. iMiller et all ()2003f ) performed their 
analysis for five Kerr BHs with different angular momenta, 
and for each case they binned the Boyer-Lindquist time 
component of the energy-momentum deposition rate 
4-vector reaching r ~ 50M into 20 polar angle bins. The 
resulting energy-momentum deposition rate per unit solid 
angle (as a function of polar angle) peaks along the surface 
of a cone centered on the rotation axis. The cone has a 
half-opening angle of 7r/4, because the spatial components 
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of the energy-momentum deposition rate 4- vector are tilted 
by about 45 degrees towards the rotation axis near the 
surface of the disk. The total (integrated over polar angle) 
energy-momentum deposition rate approximately varies 
linearly with the spin of the BH, the energy deposition 
rate of a maximally rotating Kerr BH being about twice 
that of a non-rotating Schwa rzschild BH. This i s in rough 
agreement with the results of IJaroszvnskil (jl993l ). 

The previous studies discussed above were con- 
cerned with a restricted set of specific, idealized cases : 
spherical neutrinospheres (jSalmonson fc Wilsonl 119991 : 
lAsano fc Fukuvamal |2000| ). geo metrically infinitely thin 
disks and hence neutrinos pheres (jAsano fc Fukuvam'al2000l 
I2OOII: iMiller et all I2003D . and specifically designed se- 
quences of torus models without the possibility to identify 
the separate influence of different com ponents and prop- 
erties of the neutrino emitting system (jJaroszvnskilll993l 
Il996l) . In the following we present a more comprehensive 
and systematic study of energy-momentum deposition by 
z/iz-annihilation. We analyze how the neutrino distribu- 
tion and the resulting i^£/-annihilation arc influenced by 
GR effects, the geometry and the properties of the neu- 
trinosphere, and the mass and spin of the central BH. To 
this end, we examine the most likely configurations oc- 
curring in compact astrophysical systems: spherical neu- 
trinospheres present in NSs, and thin disks and thick tori 
surrounding stellar-mass BHs. We explicitly remark that 
the absolute values of the neutrino luminosity and of the 
energy-momentum deposition rate by i/j/-annihilation do 
not play an important role for the discussion in this pa- 
per. Due to the simplicity of the neutrino source models 
considered here, our numbers for these quantities should 
not be interpreted quantitatively. They sensitively depend 
on the location and temperature of the neutrinospheres, 
which must be expected to be different in detailed hydrody- 
namic simulations that include some treatment of the neu- 
trino transport. However, for the purpose of comparing the 
different effects, on which this paper focuses, only the rela- 
tive changes of the neutrino luminosity and z^i/-annihilation 
matter. 

The paper is organized as follows: In Sect. [2] we dis- 
cuss the theoretical fundamentals for calculating the vv- 
annihilation rate in a given Kerr spacetime, and how we 
construct equilibrium accretion tori surrounding Kerr BHs. 
In Sect. [3] we give a description of the numerical implemen- 
tation of the ray-tracing algorithm used to calculate the 
annihilation rate. The results of our parameter study are 
discussed in Sect.|4l and the conclusions of our work are pre- 
sented in Sect.O Finally, technical details concerning the 
calculation of the annihilation rate are given in App. j^ and 
convergence tests performed on our ray-tracing algorithm 
are presented in App.lBl 



2. Theoretical fundamentals 

Unless stated otherwise, we use geometrized units through- 
out this paper, so that c = G = 1, where c is the speed of 
light in vacuum and G is Newton's gravitational constant. 
Greek and Roman indices denote spacetime components (0- 
3) and spatial components (1-3) of 4-vectors, respectively. 
The signature of the metric is chosen to be (-1-,— ,—,—), 
and 3-vectors arc denoted by a bar above the respective 
symbol, e.g. A. 



z 

— X 




Fig. 1. Ray-tracing of neutrinos in a Kerr BH spacetime. At 
every point x where the annihilation rate is to be calculated, 
geodesies (red curves) arriving from random directions are 
traced back until they hit the neutrinospherc (blue torus). 
The computational grid is marked by the small black dots 
and its boundary by the blue circle. See the electronic edi- 
tion for a color version of the figure. 

2.1. Calculation of the annihilation rate 

In order to compute the deposition of 4-momentum 
Qf (t, x) by the annihilation of neutrinos and antineutri- 
nos of flavor i with i g {e, fi, r} into e^e^-pairs at a spatial 
point X per unit of time and unit of volume, we follow a 
formalism very close to that of iMiller et all (j2003D . 

In flat spacetime the local annihilation rate is computed 
from the Lorcntz invariant neutrino and antineutrino phase 
space distribution functions fi,^ = f^.{t,x,p) and /p. = 
/p. {t,x,p') (for their exact definition see AppendixE]) by 
the following integral 

Qt - j dVVA" {P.P') MM , (1) 

where the Af are functions defined in Appendix [X] They 
are base d on a covariant gen eralization of the expressions 
given m iRuffert et all (jl997[ ). Note that although we are 
concerned only with time-independent models here, we in- 
clude the argument t in the distribution functions f,^. and 
for the sake of generality. In the following we omit the 
indices of /. 

Except for minor r edefinitions Eg. (fTj) can also be used 
in curved spacetimes (jMiller et al.ll2003l ). because the an- 
nihilation process, which is a microphysical phenomenon, 
happens so rapidly and on such small length scales that the 
effects of stellar gravitational fields can be safely neglected. 
Gravity affects only the propagation of the (anti)neutrinos 
between their emission (or emergence from the neutri- 
nospherc) and their annihilation. 

In our models the gravitational field, which leads to ray 
bending and redshift, is provided by the central Kerr BH of 
mass M and (dimensionless) angular momentum parameter 
a = J/M^ (where J is the angular momentum of the BH, 
and < a < 1), whose metric is given in Boyer-Lindquist 
coordinates {t, r, 6, (f>) by 

ds2 = dt^ + d(/)2 + 2gt^ dtdcp + grr dr^ + gee dO"^ (2) 
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with 

gtt 

94"P = 
9t<j> = 

9rr 

gee = 
where 



1- 



2Mr 



7 

A ' 

V, 



sm 



{aMf 



A = - 2Mt 

^ + (aMf cos^ 9 . 

Let us now consider an observer located at the point 
X where the annihilation happens (Fig.[T]). We will refer to 
this observer as the local observer, and the quantities mea- 
sured in his frame are denoted by the subscript 'L'. The lo- 
cal frame is defined by an orthonormal base {et, e^, eg, e^} 
such that a local observer is at rest in the global (r, 9, (p) co- 
ordinate system, i.e. its 4- velocity fulfills = (the result- 
ing orthonormal base vectors are given in App.[X|. Hence, 
one cannot calculate the annihilation rate inside the ergo- 
sphere, where no observers can be at rest. This restriction 
could be lifted by using observers dragged along by the BH, 



i.e. for observers with ' 



9al3U' 



0. However, the energy 



released inside the ergosphere will end up in the BH, and 
thus it is of no relevance for a GRB. We therefore exclude 
the region inside of the ergosphere from our analysis. 

In the immediate surroundings of a point x a curved 
spacetime can be approximated by a tangential flat space- 
time. Thus, a local observer can use Eq. ([1]) to calculate the 
annihilation rate. Consequently, one has to generalize only 
the momenta p in Eq. ((T|) to pl, in order to compute the 
annihilation rate Qf in the frame of a local observer. 

For doing so, one also needs to calculate the 
(anti)neutrino distribution functions in the frame of the lo- 
cal observer / (f, We take first the 4-momentum vec- 
tor components pf^ = (i?L,PL), wittQ ^-l = measured 
in the local frame and transform them into the components 



measured in the global Boyer-Lindquist frame. Hence, the 
spatial components p = p°- are functions depending on pl, 
such that we get a new distribution function f{t,x,p) = 
f (t, x,pi, (p)). Assuming now that the (anti)neutrinos do 
not experience collisions, their propagation is described by 
the collisi onless Boltzmann equation in curved spacetime 
(see, e.g.. lMisner et al.lll973l) 



d/ (f(A),x(A),p(A)) 
dA 



= 0, 



(3) 



where A is a parameter along the neutrino path. This equa- 
tion implies that the distribution function is constant along 
particle trajectories. Since we consider (anti)ncutrinos to be 



In the local frame 



3\2 



massless, their propagation takes place along null geodesies. 
The latter are given by 

dx^ 
dA 
dp" 

"dA 



P 



J \ — ^ l3fP P 



with gapp^p^ = (see, e.g.. iRiffert et al.lll998D . To com- 
pute the distribution function / (t, x,p) of either a neutrino 
or an antineutrino at the spacetime location {t, x) of mo- 
mentum p, we simply use the above equations to trace the 
path of that (anti)neutrino back in time until it hits its 
neutrinosphere (Fig.[T]). In most cases this will never hap- 
pen, and then f(t,x,p) = 0. In the other cases, we get 
an emission event at {tE,x-E) with a neutrino 4-momentum 
= {Ee,Pe)- Hence, using Eq. ^ this implies 



f{t,X,p) = f{tE,XE,PE) 



(4) 



Note that / {tE, xetPe) is the distribution function at the 
emission point in the global frame, and it has to be re- 
lated to the corresponding function in the frame comov- 
ing with the neutrinosphere. This function has the form 
of a black body for fermions (the chemical potential is ne- 
glected, which is justified for the typical conditions met in 
the systems we are interested in), and therefore it depends 
only on the comoving frame energy Eq. Taking into ac- 
count that the 4- velocity uq (t, xe) of the neutrinosphere is 
equal to the time-like base vector of the comoving frame, 
the comoving frame energy Ec is the projection of the 4- 
momentum pE onto uc (t, xe)- 

Eq ^ Pe ■ uq ■ 

Hence, the rhs of Eq. ([4]) can be expressed in terms of the 
comoving frame distribution function as 



fitE,XE,PE) 



1 



1 -I- exp 



\kBTc J 



\pl\' = {piy + {ply + (pI) 



where fee = 1.381T0 ^^ergK ^ is the Boltzmann constant, 
and Tc the temperature of the neutrinosphere. 

2.2. Neutrinospheres 

The conditions at the neutrinosphere are of importance 
when calculating the neutrino and antineutrino luminosities 
and z/P-annihilation, because the neutrino emission proper- 
ties sensitively depend on the neutrinosphere temperature 
Tc, and on the radiating area. We therefore investigate var- 
ious neutrinosphere geometries, constructed in two differ- 
ent ways. In the first approach, we prescribe the location 
and the shape of the neutrinosphere for several idealized 
BH-accretion disk systems and study their influence on the 
annihilation rate. In the second, more elabora te approach 
we proceed similar to I Jaroszvnskil ()1993l Il996f ). and model 
the accretion disk as a non-selfgravitating, stationary, geo- 
metrically thick equilibrium torus rotating around a Kerr 
BH with mass M and rotation parameter a. 

The distribution of the specific enthalpy of the equi- 
librium tori is comput ed according to the method of 
lAbramowicz et all ()1978l ) for disks with uniform specific an- 
gular momentum (/). Our equilibrium torus models are con- 
structed by fixing their mass (mtor) and the value of the in- 
ner equatorial radius of the torus (see, e.g.. lFont fc Daign^ 
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[200l) . which we choose to be Tin = 4.1 A/ and whose 
specification replaces the specification of /. Once the spe- 
cific enthalpy distribution is found, we obtain the rest 
mass density (p) distribution using the assumption that 
the equation of state is barotropic. We further assume that 
the dimensionless entropy per baryon carried by photons, 
s-y = 4a^(fcBT)'^/(3p/mii). is given and uniform, tooQ. This 
assumption translates into a relation between the density 
and the temperature (p cx T^), which allows us to com- 
pute the latter. The resulting tori are radiation dominated 
for values of s-y > 1. In Sect. 14. 21 we take = 1 in the 
tori, which corresponds to total entropies per baryon in the 
range of 5 to 12 and agrees with the values resulting from 
detailed hydrodyna mic Newtonian simulations of mergers 
of co mpact binaries (jRuffert fc Janka|[l999t ISetiawan et al.l 
^QQB)- We further assume that matter in the torus con- 
sists of free protons, neutrons, e~, e"*", and photons. For 
the typical conditions in the torus, it is justified to ne- 
glect the effects of strong interactions on the equation of 
state. With the baryon rest-mass density p, electron frac- 
tion Ye = (rio- —nc+)/nB (where ne, and are the num- 
ber densities of baryons, and of electrons and positrons, 
respectively) and given, the thermodynamic state is un- 
ambiguously defined. The specification of is needed to 
determine the composition of the torus gas such that the 
neutrino opacities can be calculated, and its value is chosen 
to be = 0.1 in our paper. Therefore, an equilibrium torus 
is fully specified by the six parameters M, a, mtor, S7, r-in, 
and Fc- 

For the equilibrium torus models the neutrino an d 
antineutrino opacities are computed as in [jmAal (l200lh . 
The neutrino spectrum is assumed to be that of a black 
body for fermions with zero chemical potential, and the 
neutrinospheres are determined as surfaces surrounding 
the BH, at which the optical depth for equilibration of 
(anti)neutrinos is r = 2/3. The integration to calculate r 
is done along rays parallel to the symmetry axis, assuming 
neutrinos to have an average energy according to a ther- 
mal spectrum with the local temperature. We point out 
that this approach to estimate the ne utrinospheric surfaces 
differs from that of lJaroszvnskil (jl996D . who integrated neu- 
trino trajectories back to the points where r = 1. Our 
method nevertheless leads to ne utrinosphere geom etries, 
which are very similar to those of lJaroszvnskil (|1996[ ). 

However, the neutrinosphere concept is only an approx- 
imation. Given some matter distribution, e.g. in an accre- 
tion torus, the neutrinospheres roughly separate optically 
thin from optically thick regions. Strictly speaking, there 
exists no well defined surface where neutrinos decouple in- 
stantaneously from the background, because the neutrino 
opacities depend on the neutrino energy. It is also a crude 
simplification to assume that there exists a surface, exte- 
rior to which neutrinos stream freely, and interior to which 
they are in equilibrium with matter and diffuse out suffer- 
ing multiple scattering. Therefore our definition of the neu- 
trinospheres appears sufficiently good to discuss the basic 
features of I'P-annihilation in the vicinity of accreting tori 
in rotational equilibrium around BHs, although it is very 
approximative and does not take into account the local con- 
ditions along the curved geodesical paths of neutrinos. 



^ The constants used in the definition of the photon entropy 
are = 8n^/{15h^c^) = 2.082-10"^ erg~^ cm~^ rUu = 1.661- 
10"^* g, and h = 6.626 ■ 10"^^ ergs. 



3. Numerical implementation 

The z/i/-annihilation rate is computed on a spherical grid 
consisting of points in radial and Ng points in polar 
(6) direction, respectively. Typically, we use Ng = 100 
uniformly distributed points excluding a small conical re- 
gion around the rotation axis of the BH for numerical rea- 
sons. The radial grid is non-equidistant, stretching from 
the BH-horizon to a radius slightly larger than the most 
distant radial point of the neutrinosphere (Figs.|2l |3l and 
in]). The radial grid points are distributed according to 
r,;+i = (1 + A9) ri, where A9 is the angular grid spacing, 
and i = 1, . . . ,Nr (typically = 100). We assume axial 
and equatorial symmetry. 

For each of our models we trace A'^rays = 20000 rays back 
in time in random directions, starting from every grid point 
of the spherical grid exterior to the neutrinospheres. The 
procedure is performed independently for neutrinos and an- 
tincutrinos (except when the neutrinospheres are identical) . 
The time integration is performed with a fourth-order adap- 
tive step-size Runge-Kutta method until the neutrino tra- 
jectory hits the neutrinosphere. To determine when a neu- 
trinosphere is hit, we have developed the following mesh 
refinement algorithm. On a coarse grid of 150 x 150 zones 
we label the zones belonging to the optically thick region 
(interior to the neutrinosphere) , and we flag the zones which 
are crossed by the neutrinosphere (Fig.|T|). The latter zones 
are further refined, such that each of them is covered by 
150 X 150 fine zones. Again, we label each of the fine zones 
belonging to the interior of the neutrinosphere. The whole 
process has to be performed only once (before calculating 
the annihilation rate on the spherical grid). This algorithm 
allows one at each step of the ray-tracing to very rapidly 
check whether the optically thick region has been reached. 
In that case the last twenty steps of the ray-tracing are re- 
done with improved accuracy. The described method can 
deal with arbitrary (axisymmetric) neutrinosphere geome- 
tries, but it is very inefficient when applied to cases in- 
volving infinitely thin disks. For these cases our algorithm 
branches to a different 'hit-detection' procedure, and checks 
whether the (x, j/)-plane (i.e., the equatorial plane), con- 
taining the neutrinosphere and the thin disk, is intersected 
between the inner and outer radial edge of the disk. Most 
ray-traced neutrinos never hit the neutrinosphere, because 
either they are trapped by the BH or because they escape 
to infinity. 

4. Simulation results 

Though we included all three (anti)neutrino flavors in the 
previous theoretical discussion, we consider the annihila- 
tion of only j/o and v^. in our simulations, since the en- 
ergy deposition is dominated by the latter process. The 
reasons for that are the reduced production and emission 
of heavy-lepton neutri nos by GRB accretion tori (see, e.g., 
iRuffert fc Jankalll999t ) and the fact that the weak coupling 
constants for i^z?-annihilation of muon and tau neutrinos 
are roughly a factor of five lower than those of electron- 
type neutrinos. Thus, in the following we use the shorthand 
notation Q" for the annihilation rate measured in the 
frame of a local observer. 

In order to quantify the total amount of energy released 
by j/i/-annihilation per unit of time, we consider four energy 
deposition rates, which are all based on volume integrals of 



R. Birkl etal.: Neutrino pair annihilation near accreting black holes 



6 

the energy component Q* of the annihilation rate 4-vector. 
The first two rates are integrals of the rates measured in 
the local frames: 

Elt - / drd0d0V-det {gab)Q\ (5) 

Ell = I drdfld^V-det (gafc) Q* • (6) 

Here Vtot and V-^-p are the volume of the whole computa- 
tional grid and the volume of that part of the computational 
grid where > 0, respectively. In the 'up' region the ra- 
dial component of the momentum vector of the e~''e~-pair 
created by j/i/-annihilation directs outward. Ignoring hydro- 
dynamic effects of the produced e"''e~-photon plasma (i.e., 
adopting a 'free-particle picture'), the energy deposited in 
this 'up' region is likely to eventually reach a distant ob- 
server, i.e. it is the crucial region for our study. In contrast, 
the energy released between the event horizon and Vup is 
probably swallowed by the BH. 

We further consider the corresponding energy dcpo- 
sition rates meas ured by an infinitely distant observer 
(|Jaroszvnskilll993[ ): 



Elf^ = / drd0d0 J-det(5„^)Q*, (7) 



= / drd0d(/) J-det(g„;3)g*. (8) 

Note that while in Eqs. (O and ^ the integrands contain 
the determinant of the 3- metric (Latin indices), Eqs. ^ and 
([S]) involve the determinant of the 4-metric (Greek indices). 
The latter two energy deposition rates are generally smaller 
than the corresponding rate computed as the integral of Q* 
in the local frames, because they include the redshift due to 
the gravitational field. If GRBs are fueled by the process of 
i/;?-annihilation, Ell'°° (Eq.[51) represents a rough estimate 
for the maximum luminosity of the GRB event. 

Two other quantities relevant for analyzing our models 
are the 'local luminosity' radiated from the neutrinosphere, 

^'' = 1^^3(0)/ daT^, (9) 

' ^ J jy— sphere 

where we assume blackbody emission from an isotropically 
radiating surface, and the corresponding 'luminosity' for an 
observer at infinity, 

LT = ^Fs{0)f daguT^. (10) 

J sphere 

Here ^3(0) = dx [x^ /{e^ + 1)] da is the appropri- 
ate GR surface element of the neutrinosphere. The factor 
Qtt in Eq. (|10p accounts for the gravitational redshift. Note 
that what we call 'luminosity' here is actually the total en- 
ergy emission rate from the neutrinosphere, which is not 
the quantity directly measured by an observer at a fixed 
position, because we ignore the direction-dependent varia- 
tion of the surface area of the neutrinosphere visible to an 
observer. 

The evaluation of the quantities E'*"*'"" and Ell'°° (see 
Eqs. ([7]) and ©) requires two steps. In the first step Q* is 
calculated at the annihilation point as described in Sect. 
12.11 including aberration effects due to the neutrinosphere 



rotation (for models with 1^0). Using the obtained Q*, 
Eqs. ([7]) and ^ arc evaluated in the second step. Note that 
these two equations involve an approximation, because the 
spatial components of the annihilation rate 4-vector do 
not enter. This approximation means that the momentum 
of the i^P-annihilation plasma at the annihilation point is 
ignored in the transformation of the energy to infinity. In 
contrast to the evaluation of Q*, neutrinosphere aberration 
effects are neglected in Eq. (flOl) . because this enormously 
simplifies the evaluation and can be justified by the use of 
an approximation in Eqs. ([7]) and ^ as explained above. 

4.1. Idealized models 

In the following we investigate several idealized BH- 
accretion disk systems. For each of these models the cen- 
tral black hole of mass M and angular momentum a 
is surrounded by a constructed neutrinosphere (see also 
Sect. l2.2p . which is characterized by the geometry (thin 
disk, torus, or sphere), the position (radii), the rotation, 
and the temperature (see Table [1]). The rotation of the 
neutrinosphere is given by the uniform Lagrangian angu- 
lar momentum I = —u^/ut (jAbramowicz etall[l97l . This 
is the GR generalization of the specific angular momen- 
tum I = X X V, where x is the position and v the velocity 
In our models the prescribed geometry and position of the 
neutrinosphere limit physically reasonable values for the 
Lagrangian angular momentum I to several times the mass 
of the BH (cgs units are obtained by applying the conver- 
sion factor (Mgr/Mq) • 4.42 • 10-^^ cm^ s~^), and values of 
I K, AM approximately describe Keplerian rotation. The 
neutrinosphere is assumed to be isothermal in the comov- 
ing frame, i.e. it has a uniform comoving frame tempera- 
ture Tq. In each idealized model the four parameters of the 
neutrinosphere presented here are also used for the antineu- 
trinosphere. 

The idealized models are constructed for comparison 
with two reference models, which arc the disk models D and 
REF (see Tabled]). Both models have the same temperature 
Tc, and the disk is assumed to be non-rotating (/ = 0) for 
computing the z/i/-annihilation effects, but the masses M 
of their non-rotating black holes and their outer disk radii 
differ. For each reference model there is a series of models 
(models DN to TL4N for model D and models A2B to L5 for 
model REF, see Tabled]), which are obtained by changing 
the properties of the corresponding reference model. We 
explore the influence of the geometry (toroidal models T, 
A2B, AB, A.5B, and TEMP, and spherical models S, and 
SMS), of GR effects by studying 'Newtonian' models (disk 
model DN, torus model TN, and sphere model SN), of the 
inner disk radius (models RI4.5, and RI3), of the black hole 
angular momentum (models DAI, A. 5, and Al), and of 
the disk rotation (models DL6.5, DL4, TL4, DL4N, TL4N, 
L2.5, L4, and L5). 

For most of the idealized models the effects due to 
the rotation of the neutrino emitting surface are ignored. 
Such models correspond to the cases which have / = 0. 
Although in astrophysical systems disks and tori always 
rotate, we nevertheless consider idealized models with no 
neutrinosphere rotation here, because we want to selectively 
study the effects of disk or torus rotation on the neutrino- 
antineutrino annihilation rate. Later (see Sects. 14.1.31 and 
14.1. 4p . we will demonstrate that the rotation of the neu- 
trino emitting source has a significant effect only on the 
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Table 1. Some properties of the simulated idealized models, whose isothermal neutrinosphere and antineutrinosphere 
are assumed to coincide and to have the same temperature. The quantities given in the columns of the table from 
left to right are: the mass M and the (dimensionless) angular momentum parameter a of the BH, the geometry of 
the (anti)neutrinosphere, the inner and outer radius (disk/torus) or just the radius (sphere) of the neutrinosphere; for 
toroidal neutrinospheres the ratio of the radial to the vertical diameter of its elliptical meridional cross section is given, 
too. The other quantities are the Lagrangian angular momentum I, measured in terms of the BH mass Mgr of the GR 
model, i.e. in case of a GR model we have Mgr = M ^ 0, whereas in case of a Newtonian model [M = 0) the value 
of the Lagrangian angular momentum of the corresponding GR model is used (cgs values of I are obtained by applying 
the conversion factor (Mgr/M©) • 4.42 • lO""^^ cm^ s~^), the comoving frame temperature Tc of the neutrinosphere (in 
the special model TEMP there is no emission outside of a cylinder around the system axis with radius 8M), the total 
neutrinosphere luminosity L'^ ~ + for an observer at infinity, the total energy deposition rate, , as the 
integral over the total volume of the local energy deposition rate Q*, the total energy deposition rate i?"? in the 'up' 
region (defined by > 0), these total energy deposition rates measured by an infinitely distant observer, and 
and the efficiencies q°^"°° = and = respectively. 
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spatial distribution of the annihilation rate but almost no 
influence on volume integrated results. 

The annihilation rate 4- vector of the series of models 
D to TL4 is shown in Figs. [2]and[3l Si milar plots can also be 
found in Fig. 5 of lMiller et al.l (|2003f ). and there especially 
the upper left panel (a = 0.0, 'MDR') nearly looks like 
our plot for model DL6.5 (see Fig.[3|), where the accretion 
disk rotates with the super-Keplerian value / = 6.5. The 
agreement is even somewhat better in case of model L5 (not 
visualized in our work), because in contrast to model DL6.5 
it has the same inner (6 M) and oute r (10 M ) disk radii as 
the mentioned model of lMiller et al.l (l2003l). However, we 
were not able to reproduce the models of lMiller "eran(l2003l) 
exactly, because of the lack of information about the precise 
parameters of their models, particularly the rotation law. 
We als o point out that in the pictorial representation of 
Fig. 5 of lMiller et al.l (j2003[) . the local value of the azimuthal 
integral of the 'MDR' is provided, while our figures show 
the local value of Q'. If instead of Q* we also compute the 



azimuthal integral of Q*, the contours of, e.g., model DL6.5 
be come vertica l in th e vicinity of the rotation axis just like 
in iMiller et al.l (l200l . 

4.1.1. Influence of general relativity 

In order to study the influence of GR effects, we take the 
thin disk model D, the torus model T, and the sphere 
model S, and compare these models with their correspond- 
ing Newtonian cases DN, TN, and SN (see Table [1]). The 
Newtonian models have the same neutrinosphere as the cor- 
responding GR models, however, we assume M = 0, i.e. we 
neglect the influence of gravity on the neutrino propagation 
and transformation between observers. The computational 
grid for a Newtonian model is identical to that of the corre- 
sponding GR model, i.e. in particular, the inner radial grid 
boundary is located at the same radius in both models (just 
above the horizon of the GR model). All three Newtonian 
models DN, TN, and SN are assumed to have the same 
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Fig. 2. Annihilation rate 4- vector measured in the frame of a local observer for models with idealized neutrinospheres 
(see Tabled]). The location of the neutrinosphere, which coincides with the antincutrinosphere, is marked by a dashed 
line. The value of the energy component Q* is color-coded, while the spatial vector Q is visualized by showing the spatial 
velocity vector v = Q/Q* , whose component coplanar (perpendicular) to the displayed x — z-plane is given by the black 
arrows (filled circle). The length (area) of the arrow (filled circle) is a linear measure of the size of the component, the 
maximum value (c = 1) being represented by the arrow (filled circle) in the right upper corner of each panel. Note that 
for the neutrino evaluation the neutrinosphere is assumed to be non-rotating. Therefore the perpendicular component of 
V is zero here. The big black circular region in the left lower corner of the left two panels represents the BH. In the right 
two panels the corresponding region is chosen to be white, because these panels show Newtonian models where the GR 
effects of the BH were disregarded. Finally, in every panel a solid white line separates regions with positive from those 
with negative . Note that in model DN the white region along the equatorial plane (x-axis) is caused by the saturation 
of the color scale due to a very low energy deposition rate. See the electronic edition for a color version of the figure. 



surface area and temperature, and therefore also have the 
same total (neutrino plus antineutrino) luminosity , i.e. 
iiy,DN = iiy,TN = -^ly.SN = 3.9 • lO^^ergs"^ (see Table [T]). 
Although the GR models D, T, and S have the same met- 
ric (M = 2MQ,a = 0), a local comoving observer does 
not measure the same surface areas, because the surfaces 
of these three models do not coincide such that they are dif- 



ferently affected by the curvature of spacetime. Therefore, 
despite having the same temperature {Tc = 5 • 10^° K), 
models D, T, and S have a different local luminosity (see 
Eq.[5]). Note that in the Newtonian models the energy emis- 
sion rate from the neutrinosphere and the energy deposition 
rates do not depend on the observer, because there is no 
gravitational redshift effect. 
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Fig. 3. Same as Fig. [51 but for a different set of models with idealized neutrinosphere geometries (see Table [T]). Note 
that model DAl contains a central rotating BH (with a = 1) and therefore possesses an ergosphere outside of the event 
horizon. This ergosphere is marked by the gray area. See the electronic edition for a color version of the figure. 




Fig. 4. Energy component of the annihilation rate 4-vector along the z-axis for models with idealized neutrinosphere 
geometries (left panel; see Table [T]), and for models where the neutrinospheres are calculated from equilibrium torus 
models (right panel, see Table [3]). The values are normalized to the energy component of the annihilation rate 4-vector 
along the z-axis of models D and E.01.05, respectively. 
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Fig. 5. Energy deposition rate E^p'^'°° on a sphere of radius 200 M per unit of Boycr-Lindquist polar angle 6, obtained 
by taking into account the energy transport to this radius as described by the annihilation rate 4-vector field Q". The 
left panel shows the results for idealized disk models with different values of the Lagrangian angular momentum /, and 
the right panel the results for the first four equilibrium torus models of Table |3l 



In the local frame GR effects increase the energy de- 
position rate E^p of models D, T, and S with respect to 
their Newtonian counterparts (see Table [J) • The increase 
is smaller when restricting the comparison to the 'up' re- 
gion (where the radial component of the momentum depo- 
sition vector is positive, > 0). But even in the 'worst' 



case (toroidal model T) GR effects enhance E'"? by about 
40%. The optimal geometry to release energy in the sur- 
roundings of the BH is that of a thin disk, where the GR 
enhancement of the energy deposition rate in the 'up' re- 
gion is a;"p = 2.3 for a local observer and where the en- 
ergy deposition rate in the 'up' volume is larger than for 
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Table 2. Comparison of idealized models: For each pair 
of models x^p = E,^[, (Ml) /E^[^ (M2) denotes the ratio of 
the total annihilation rates in models Ml and M2. The 
superscript 'tot' indicates that this ratio is obtained by in- 
tegrating the local energy deposition rates over the whole 
grid, while for the superscript 'up' the integration is re- 
stricted to the 'up' region, where the energy released by 
i^i/-annihilation is not trapped by the BH and may eventu- 
ally reach a distant observer. The ratio is computed in the 
local frame or in the frame of an infinitely distant observer. 
The latter case is indicated by the additional superscript 
'oo'. 
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toroidal and spherical geometries. These results also hold 
for a distant observer (Tabled]). The energy deposition by 
z/p-annihilation is therefore enhanced in GR models com- 
pared to the Newtonian treatment, despite the fact that the 
luminosity LJf for an observer at infinity is smaller (Table 
[1]) . This in the first moment counterintuitive behaviour can 
be understood by the relativistic effects, which account for 
a more intense neutrino radiation field close to the emitting 
surface, i.e. the local luminosity, which is not reduced by 
gravitational redshift, is always higher or at least equal to 
the Newtonian counterpart of a GR model. 

The enhancement of the energy deposition rate due to 
GR effects depends strongly on the position, and becomes 
negligible at sufficiently large distances from the BH (Fig. 51 
left panel), where relativistic effects are unimportant. In 
case of the disk models (D and DN). Fig. [5] (upper panels) 
and Fig.d] (left panel) show that GR effects enhance the 
energy deposition rate near the symmetry axis {0 < 30°) 
by a factor of > 10 up to r ^ 10 km, and they still increase 
the energy release by a factor of ~ 2 at r ~ 30 km. 

4.1.2. Influence of the neutrinosphere geometry 

According to Table [2l the energy deposition rate by lyp- 
annihilation in the 'up' region measured by an observer at 
infinity is largest for a neutrinosphere having the geometry 
of a thin disk (model D; = 2.3 • lO'^^erg s'^; Table 

[1]), while it is smallest for a neutrinosphere of spherical 
shape (model S; E]^?'°° = 0.83 • lO^^ gj.gg-i. Table|T]). This 
finding holds for the corresponding Newtonian models (DN 
and SN), too. This is a remarkable result, as in the spherical 
models Q'' > everywhere on the computational grid (for 
models S and SN the ^/I'-annihilation integrals in the 'up' 



regions are equal to the corresponding integrals in the 'tot' 
region, respectively; see Table [T]). 

The influence of the neutrinosphere geometry on the 
energy deposition rate can also be studied by considering 
toroidal models with different meridional cross sections. All 
our toroidal models have an elliptical meridional cross sec- 
tion characterized by its eccentricity (Tabled]). Model REF 
(an infinitely thin disk model with an eccentricity e = 1) 
serves as the reference model, which is compared to models 
A2B (e = \/3/2, oblate), AB (e = 0, circular) and A.5B 
(e = \/3/2, prolate), respectively. These four models have 
the same Newtonian luminosity , and the center of their 
meridional cross sections or the mid-points of their disks are 
all located at Vc ~ 18 km in the equatorial plane. Compared 
to the thin disk model REF, the other torus models are all 
less efficient in depositing energy in the 'up' region (Table 
12]) ■ By inspection of Tabled] one recognizes that the energy 
deposition rate in the 'up' region measured locally and by 
a distant observer E^^'°° decreases when the meridional 
cross section becomes more prolate. It is smallest for model 
A.5B, which is the most prolate model. At the same time 
this model is the most efficient one in the total energy re- 
lease (2;*p*'a°5b = > !)■ This at first glance surprising 
behavior can be explained by the fact that a thin disk emits 
a large fraction of the neutrinos in directions parallel to the 
symmetry axis, and hence just above and below the thin 
disk. Consequently, i^i^-encounters are less frequent close 
to the BH than in a source of toroidal shape. Although 
infinitely thin disks are a mathematical idealization more 
than a natural case, our results nevertheless imply that a 
distant observer receives more energy from systems which 
have a more oblate structure of a toroidal accretion disk. 

Another geometrical aspect of relevance is the size of 
the innermost radius of the disk (rin). To study the depen- 
dence of the energy deposition rate on this parameter, we 
consider the thin disk models RI4.5 (nn = 4.5A/) and RI3 
(rin = 3M), which both have an innermost radius smaller 
than that of model REF (nn = 6A/), as can be seen from 
Table [1] but have the same Newtonian luminosity and 
Newtonian surface area. Independent of the region of in- 
tegration ('up' or 'total') and of the observer (local or dis- 
tant), we find that the smaller is Hn, the larger is the energy 
released by i^i/- annihilation (Table |2]). 

Finally, wc consider model TEMP (Tabled]) to further 
investigate the influence of the neutrinosphere geometry 
on the energy deposition rate. This model is identical to 
model AB, except that in model TEMP only the inner 
half of the torus (up to a distance of 8M from the sym- 
metry axis) is assumed to emit neutrinos. Thus, in order to 
have the same Newtonian luminosity in both models, 
the neutrinosphere temperature has to be increased from 
Tc = 5- 10^° K to Tc = 6.1 • 10^° K, because the emitting 
surface is smaller. The more concentrated neutrino radia- 
tion field therefore leads to a strong overall increase of the 
energy deposition rate (xi/p > 3; Table [2]), in analogy to 
the models of the previous paragraph, where the innermost 
disk radius was reduced. 

4.1.3. Influence of the angular momentum of the BH 

The dependence of the energy deposition rate on the an- 
gular momentum of the central BH is studied using model 
DAl, which consists of a disk surrounding a maximally ro- 
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Table 3. Properties of equilibrium torus models, where the location of the t'o-sphere differs from that of the Po-sphere. 
Besides the quantities already defined in Table [1] the models here are characterized by the mass mtor of the accretion 
torus, and the mass Aftor and angular momentum ator of the BH used to calculate the structure of the accretion torus. In 
most models the values of the last two quantities are equal to the values of the mass M and angular momentum a of the BH 
that we use to determine the neutrino geodesies. However, in the two models E.Ol.ld and El. Id we have chosen otor 7^ a, 
in order to discriminate between the influence of the black hole rotation on the accretion torus and on the neutrino 
trajectories (see Sect. 14. 2p . For the Newtonian model E.Ol.lN the table shows M ^ Mtor (and also a ^ ator), because the 
neutrino trajectories are evaluated without GR effects. All models in this table have the same photon entropy = 1, 
the same inner equatorial torus radius r^^ = 4.1 M, and the same electron fraction Ye = 0.1 (see Sect. 12. 2p . Moreover, 
the effects of the torus angular momentum I are included in the evaluation of the i^i/-annihilation. Note that I is derived 
from the set of six quantities (M, a, mtor, S7, Hn, ^o) that fully specify a torus in rotational equilibrium. 
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tating Kerr BH (a = 1; Tabled Fig.[3l). Compared to the 
corresponding non-rotating model D (a = 0), we do not find 
an enhancement of the energy deposition rate, independent 
of the region of integration and of the locatio n of the ob- 
server (Tabled)). This result differs from that of lMiller et al.l 
(j2003t ). who obtained a factor of ^ 2 larger amount of de- 
posited energy when comparing thin accretion disks orbit- 
ing a maximally rotating and a non-rotating BH, respec- 
tively. This discrepancy arises because we set = 6 M 
both in model D and DAI, in order to just test the in- 
fluence of BH rotation. In contrast, the disks considered by 
iMiller et al.l extended from a fixed outer radius Tout = 10 A/ 
in towards the BH, i.e. to the innermost stable circular or- 
bit. This results in disks with smaller values of as a ^ 1 
than in our study, and hence, as we have shown in the last 
Section, an increased energy deposition rate is expected. 
Table [T] also contains the two models A. 5 and Al, which 
are equal to our second reference model REF (a = 0, and 
Mbh = SMq, in contrast to Mbh = 2M0 for model D), 
except that the black hole rotates with a = 0.5 and a = 1, 
respectively. We find that for every considered BH mass the 
value of depends only weakly on the value of a. 

4.1.4. Influence of the accretion disk rotation 

In all idealized models discussed up to now with the shape 
of an infinitely thin disk, a torus, or a sphere, respectively 
(Table [1]), the rotation of the neutrinosphere was neglected 
in calculating the :^P-annihilation. We now address the con- 
sequences of relativistic aberration connected with such ro- 
tation on the energy deposition rate. Our idealized rotating 
disks are assumed to have a uniform Lagrangian angular 
momentum /. In Model DL6.5 the disk is rapidly rotat- 
ing with the super-Keplerian value I = 6.5 M (correspond- 
ing to a cgs value of the specific angular momentum of 
/ = 5.7 • 10^^ cm^ s-i for the IMq BH), but otherwise the 
same parameters as in the non-rotating model D (/ = 0) 



are used. For an observer at rest, the (anti)neutrinos are 
preferentially emitted tangentially to the disk in the direc- 
tion of rotation. This behavior results in the existence of 
a non-vanishing Q'^ component and is reflected by the in- 
creased size of the filled circles in the vicinity of the disk in 
Fig. [3l Consequently, neutrinos and antineutrinos collide 
less frequently with large relative angles, and in the 'tot' 
region of model DL6.5 less energy is deposited than in the 
non-rotating disk model D {x^°^ < 1; Table [2]). However, if 
we consider the 'up' region, this effect is compensated by 
the slightly increased size of the 'up' region in model DL6.5 
~ 1), especially in the vicinity of the accretion disk 
(in model DL6.5 the white line crosses the equatorial plane 
at a smaller radius than in model D, see. Figs. O and [3]). 
The results also hold for an observer at infinity. Hence, in- 
dependent of the choice of the observer, the influence of the 
disk rotation on the integrated energy deposition rate in the 
'up' region is small, which justifies the remarks of Sect. 14.1] 
regarding the choice of non-rotating reference models. This 
is also true for intermediate values of / (see Table [1]), i.e. 
in case of the sub-Keplerian rotation of model L2.5, mod- 
els DL4 and L4, whose accretion disk approximately ro- 
tates at Keplerian speed, and the super-Keplerian model 
L5 (the three models L2.4, L4, and L5 are not based on 
model D, but on model REF). Also in case of a toroidal 
neutrinosphere instead of a disk, the effect of the neutri- 
nosphere rotation on the total annihilation rate measured 
by an observer at infinity is only small (see the values of 
for models T and TL4 in Table [J). 

We have also studied the influence of the neutrinosphere 
rotation in case of Newtonian models. For this purpose, 
we have calculated models DL4N and TL4N (Table [T]). 
When GR effects are ignored we again find that the neu- 
trinosphere rotation does not have an important impact on 
the total annihilation rate measured at infinity (compare 
for models DL4N and DN, and TL4N and TN in 
Tabled]). 
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Fig. 6. Local annihilation rates for neutrinosphere geometries calculated from four equilibrium torus models (see Table 
[3]). The visualization method is the same as in Figs. [5] and O but here the neutrinosphcres for ly^ (blue, long dashes) and 
(red, short dashes) do not coincide. Note that the accretion tori of the four models shown here rotate, however the 
spatial components of the annihilation rate 4- vector perpendicular to the displayed x — z-planc are very small and visible 
only in the close vicinity of the black hole. See the electronic edition for a color version of the figure. 



Despite the weak dependence of the integrated values of 
the energy deposition rate in the 'up' region on the specific 
angular momentum of the disk, the spatial distribution of 
the energy deposition rate changes significantly (compare 
the plots of models D, DL4, and DL6.5 in Figs.[l] and [3]). 
The faster the disk rotation is, the larger are the spatial 
components of annihilation rate 4-vector tangential to the 
disk. In addition to that, the energy deposited in the vicin- 
ity of the symmetry axis decreases, which is especially vis- 
ible in model DL6.5. 

For the non-rotating model REF and the two rotating 
models L2.5 and L5 the left panel of Fig. [5] shows the en- 
ergy dE^p''°"/d6 arriving at a sphere of radius 200 A/ per 
unit of time and per unit of Boyer-Lindquist polar angle 
9. Numerically, these results are obtained by closely fol- 



lowing the approach of iMiller et al.l (|2003f ). We consider all 
points of our numerical grid and for each grid point we 
transport the calculated annihilation rate 4-vector par- 
allel to itself, i.e. wc follow geodesies in the direction of 
the 4-vector, starting at the annihilation point. Excluding 
the fraction of parallelly transported 4- vectors that eventu- 
ally hit the black hole or the neutrinosphere, the remaining 
4-vectors Q"' are transported until they hit the sphere of 
radius 200 M (centered in the system origin). For each such 
4-vector, we take the initial energy component Q*, i.e. be- 
fore the parallel transport, and multiply it by the factor 
27rArA0^-det (gap), in analogy to Eqs. and ©. This 
way we take into account the 3-volume and gravitational 
rcdshift. Finally, depending on the 0- angle of the 4- vector 
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Q" arriving at the sphere, we create the plots of Fig.O by 
using 20 equally spaced 0-bins. 

In case of the rotating thin disk models L2. 5 and L5 (see 
Fig.[S]) we obtain results very similar to those of lMiller et alj 
()2003f ). namely the dominant contribution to dE^'^^'°° /dO 
is found to be off-axis and scaling their results to cgs units 
they quantitatively agree with ours (compare Fig. 8 in their 
work with the left panel of our Fig. [5]). However, the position 
of the peak strongly depends on the disk rotation, and for 
non-rotating models the bulk of the energy is deposited 
in the vicinity of the system axis. In addition to the disk 
models REF, L2.5, and L5, we have also considered the 
non-rotating torus model AB and rotating idealized torus 
models (not shown in this work), and we find the same 
general behaviour as in case of the disk models. 

4.1.5. Efficiency 

The efficiency of converting radiated neutrino energy into 
e"^e~-pairs by i^S'-annihilation is given by the quantity 
q^^'°° = where L"^ is the sum of the neu- 

trino and antineutrino luminosities at infinity (see Table 
[T|). Values as large as a few tenths of a percent can be 
reached for (7"^'°°. Neutrinospheres having the geometry of 
a thin disk are the most efficient ones in converting neutrino 
energy into e~'"e~-pairs, while neutrinospheres of toroidal 
shape have the smallest efficiency (see the models D, T, 
and S in Table [T|). If GR effects are not taken into account 
in calculating the energy deposition rate by i^P-annihilation 
(see models DN, TN, and SN), thin disks are still most effi- 
cient, however, now the spherically shaped neutrinosphere 
is least efficient. Higher disk luminosities increase the depo- 
sition efficiency (see the model REF and all models below in 
Table[T|). The efficiency monotonically depends on the inner 
radius r^^ of accretion disks (see models REF, RI4.5, and 
RI3 in Table [T]). When rjn is reduced, the neutrinosphere 
moves closer to the event horizon of the BH. Although, 
this way the neutrino luminosity measured by an observer 
at the locations of i^P-annihilation rises (not shown in Table 
[1]), the neutrino luminosity at infinity nearly remains 
the same due to the increased gravitational redshift. Since 
the energy deposited at infinity E]!,p'°° increases, the effi- 
ciency becomes larger for smaller values of the inner disk 
radius r-m. An analogous effect can also be seen in the torus 
model TEMP (see Table [T]) , where compared to the torus 
model AB the accretion torus and Newtonian luminosity 
are the same, but the radiating region is closer to the 
BH. In contrast, up to two significant figures the efficiency 
exhibits no dependence on the angular momentum of the 
black hole (see, e.g., model A. 5 compared to model REF) 
and of the disk (see models L2.5, L4, and L5 in Table [1]). 

4.2. Equilibrium models 

Next we discuss the i^P-annihilation rate in the vicinity of 
the neutrinospheres of equilibrium accretion tori computed 
as described in Sect. 12. 2] In all of our equilibrium torus mod- 
els we keep s-y, rjn, and Y^. fixed, and therefore we will con- 
sider only the influence of the three parameters {M, a, TOtor) 
in the following paragraphs. 

In contrast to the idealized models, the neutrinosphere 
and the antineutrinosphere do not coincide in the equilib- 
rium models (Fig.lH]), because there is a larger number den- 
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Fig. 7. Temperature profiles of the neutrinospheres of the 
models in Fig. [S]in the comoving frame versus the angle t?, 
which is defined as follows: If Tc is the arithmetic mean of 
the innermost (nn) and outermost (^out) radius where the 
neutrinosphere intersects the equatorial plane, then given 
a point C on the neutrinosphere the angle i) is enclosed by 
the legs re Tout and TcC, respectively. The blue (red) lines 
denote the temperature of the neutrinosphere (antineutri- 
nosphere) . See the electronic edition for a color version of 
the figure. 



Table 4. Comparison of the equilibrium torus models. The 
table is structured like Tabled but the models shown here 
are models listed in Table [31 
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sity of free neutrons than free protons in the equilibrium 
tori. Hence, the opacity for neutrinos and antineutrinos dif- 
fers, and consequently the locations of the neutrinosphere 
and antineutrinosphere, too. Since the matter in the equi- 
librium tori is neutron rich, the antineutrinosphere is al- 
ways completely interior to the neutrinosphere. On aver- 
age it also has a higher temperature (Figs. [6] and [7|), i.e. 
Lp > L^. Neither the neutrinosphere nor the antineutri- 
nosphere are isothermal as in case of the idealized models 
(Fig. [7]). Moreover, the torus temperature now depends on 
the BH spin and increases, the closer the torus is located 
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around the BH. This effect is ignored in the ideahzed mod- 
els. 

The energy deposited in the region between the neutri- 
nosphere and the antineutrinosphere cannot be computed 
by our approach (hence the corresponding region is white in 
Fig. [S]), because in that region transport physics matters, 
i.e. neutrinos do not propagate freely. 

In order to analyze the influence of the torus mass mtor, 
let us consider the equilibrium models (Table [3]) E.01.05, 
E.01.25 and E.01.5, all of which orbit a slowly rotating BH 
with a = 0.01, and models E.8.05 and E.8.5, which gird a 
rapidly rotating BH with a = 0.8. From Table[3]we infer (i) 
that a larger torus mass strongly increases the vv-cwcigy 
deposition rate, because of a higher neutrino luminosity, 
and (ii) that this effect is stronger for lower values of a. The 
increase of the energy deposition rate is almost linear with 
mtor for high values of a, but grows like m^'^^ for a = 0.01. 
Thus, for a distant observer the energy released in the 'up' 
region can be roughly fitted by Elf°" cx (E.01.05) • 

(mto./M0)O-4-+O-6. 

On the other hand, comparing models with the same 
torus mass but with BHs that have different a, e.g. models 
E.01.05 and E.8.05, shows that larger values of a lead to a 
larger energy deposition (x"p °° =2.6; Table [3]). The effect 
of a on the deposition rate is about a factor of 2 smaller 
for higher torus masses (compare models E.8.5 and E.01.5 
in Table SI). 

What causes the increased energy deposition rate in case 
of faster rotating BHs? There are two possibilities: the ro- 
tation rate parameter a influences (i) the geometry of the 
accretion torus and hence the neutrinosphere, and (ii) it 
affects the (anti)neutrino trajectories through frame drag- 
ging. In order to determine the importance of both possibil- 
ities, we consider artificial test models where we construct 
the torus model for a rotating BH but then disregard the 
BH rotation for computing the neutrino geodesies, or vice 
versa. This means effectively that we consider two differ- 
ent rotation parameters for the same model. The first one, 
Otor, is used to calculate the equilibrium tori, and hence 
the neutrinospheres. The second one, a, is used to calcu- 
late the GR effects on the propagation and annihilation 
of neutrinos and antineutrinos. Note that the smaller of 
both values is used to set the inner edge of the computa- 
tional grid, which is thus located slightly outside the larger 
of both horizons. Such models are actually inconsistent, be- 
cause a change in the BH spin automatically translates into 
a structural change of the accretion torus. We nevertheless 
consider them in order to be able to discriminate the influ- 
ence of the BH rotation on the !/i?-annihilation through its 
effects on the neutrino trajectories from those effects that 
result from changes of the torus properties around rotating 
or non-rotating BHs. 

For each of the two parameters, a and Otor, we con- 
sider two extreme values, which leads to the four models 
E.01.1, E.Ol.ld, El. Id, and El.l, respectively (Tabled. 
As in case of the idealized models (Sect. HTTj) . increasing 
the BH's rotation rate from a — 0.01 to a = 1 does not 
yield any appreciable increase of the energy deposition rate 
(a;"p'°° ~ 1) for both values of ator (Tabled]). However, a 
substantial increase (^ 60%) of the energy released in the 
'up' region (as seen by a distant observer) is found when 
increasing ator from 0.01 to 1. Thus, we can unambiguously 
conclude that the increase of the dimensionless angular mo- 
mentum of the BH yields a larger energy deposition by vD- 



annihilation essentially exclusively because it allows for a 
smaller innermost radius of the equilibrium accretion torus. 

Concerning the influence of the GR effects, a com- 
parison of model E.01.1 with its Newtonian counterpart 
E.Ol.lN (Table O shows that GR effects on the neutrino 
propagation increase in case of the equilibrium mod- 

els only by ^ 20%. This enhancement of the energy deposi- 
tion rate is similar to that obtained for the idealized torus 
models (see Sect. 14. ip . 

In analogy to the idealized models, we also computed 
the energy deposition rate distribution on a sphere of ra- 
dius 200 M per polar angle 9 for equilibrium models (see 
right panel of Fig. [S]). However, in contrast to the idealized 
models, the bulk of the energy is a lways deposited nea r the 
system axis. Hence, the result of iMiller et al.l (pOOS*), i.e. 
that the main contribution to the energy deposition rate 
comes from the off-axis region, is valid only for idealized 
models (see Sect. 14.1.^ . but in equilibrium torus models 
the energy deposition along the symmetry axis is clearly 
dominant. 

The efficiency (7"^'°° increases for our equilibrium mod- 
els both with the torus mass and with a, reaching values 
of 1-3 per c e nt (T able [3]) in agreement with the results of 
iJaroszvnskil (1 19961) for his m o dels w ith specific entropies 
> 8, and with Setiawan et al.l (|2004[ ) for their models with 
the largest a-viscosity. The increase with a is a structural 
effect, because for larger values of a the torus is on average 
closer to the rotation axis and hotter, i.e. the neutrino lu- 
minosity is higher. Therefore the neutrino density is higher 
near the system axis and hence a larger fraction of neutrinos 
and antineutrinos annihilate with each other. The efficien- 
cies of the equilibrium models are ~ 10 times larger than 
those of the idealized models, because both the tempera- 
tures (Fig. [7]) and the surface areas of their neutrinospheres 
are larger. 



5. Conclusions 

The main goal of this paper is to study generic properties of 
the process of energy deposition by neutrino-antineutrino 
iyv) annihilation in the vicinity of systems consisting of a 
central stellar-mass black hole (BH) and an accretion disk 
or torus surrounding it. Assuming that the thermodynamic 
conditions in the accretion flow are such that a copious 
flux of vv-^&us can be produced, we have performed a sys- 
tematic parameter study of the influence of the mass and 
dimensionless angular momentum parameter of the BH, of 
the shape and thermal properties of the neutrinosphere, 
and of the importance of different general relativistic (GR) 
effects on the amount of energy released by i^P-annihilation. 

On the one hand we considered idealized models hav- 
ing an isothermal neutrinosphere of prescribed temperature 
and geometry, and on the other hand non-selfgravitating, 
axisymmetric equilibrium tori bound to a central BH of 
given properties. In the latter models the neutrinospheres of 
neutrinos and antineutrinos do not coincide, because they 
are computed for opacities of Vf. and that differ due to 
absorptions on free n eutrons or free protons, respectively 
(similar to the work of |Jaroszvnskilll993l . ll996[ ). Using both 
sets of models we numerically calculated the annihilation 
rate energy-momentum 4-vector. For this purpose, we con- 
structed the local neutrino distribution by ray-tracing neu- 
trino trajectories in a Kerr space-time using GR geodesies. 
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Our s tudy is a generalization of the work of iMiller et al.l 
()2003f ) to axisymmetric neutrinospheres of different shapes, 
which are not necessarily spatially coincident for i^o and Vc- 
We considered three different neutrinosphere geometries 
in our set of idealized models covering cases that may be 
encountered in astrophysical systems: infinitely thin disks, 
tori, and spheres. Infinitely thin disks are mathematical 
idealizations. However, they are investigated as the lim- 
iting case of very oblate tori. For this reason they have 
been widely used in the literature. By comparing with 
corresponding Newtonian models, where the influence of 
the gravitational field on the neutrino propagation is ne- 
glected, we have found that for an observer at infinity GR 
effects enhance the energy deposition rate by a factor of 
2 in the non-rotating, thi n disk case in agreement with 
lAsano fc Fukuvamal (|200lh . In the other two cases the in- 
fluence of GR effects enhances the annihilation rate only by 
« 25%. 

Independent of whether GR effects are included, the 
energy deposition rate that may lead to the acceleration of 
relativistic outflow is largest for infinitely thin disks. For 
more prolate tori the energy deposition rate drops in com- 
parison with a disk that has the same neutrino luminos- 
ity, because a sizeable fraction of the energy deposited by 
z/p-annihilation is swallowed by the BH. Spherical neutri- 
nospheres are the least favorable geometry for a large total 
annihilation rate. Accretion disks encountered in gamma- 
ray burst (GRB) scenarios are likely to be geometrically 
thick, in which case the geometry of their neutrinospheres 
is torus-like. For such a situation neglecting GR effects in 
the evaluation of t/i?-annihilation is a good approximation. 

We also analyzed the influence of the dimensionless ro- 
ta tion paramet e r a of the BH, and we confirm the findings 
of iMiller et al.l (j2003[ ) that modifying a without changing 
the innermost disk radius (i.e., leaving the thin disk geome- 
try unchanged) has no effect on the amount of released vD- 
annihilation energy. The independence of the results of the 
BH rotation holds both for idealized models with toroidal 
neutrinospheres and for equilibrium accretion torus mod- 
els. However, in a consistent accretion disk model, the in- 
nermost radius of the disk (close to the innermost stable 
circular orbit) will shrink as a ^ I. We have found that 
a smaller innermost disk radius leads to a substantial in- 
crease of the energy deposition rate, even when the local 
Newtonian (i.e., GR effects disregarded) neutrino luminos- 
ity is the same. This can be understood by the higher neu- 
trino density in t he vicinity o f a mo re compact disk and 
was also found bv IMiller et all (|2003D . 

Depending on the mass of the accretion torus, mod- 
els containing a maximally rotating BH (a = I) can release 
roughly a factor of 2 more energy (as measured by a distant 
observer) by :^P-annihilation than models involving non- 
rotating BHs. For a small accretion disk mass (^ O.O5M0) 
the energy release in case of a rotating (a = 0.8) BH is more 
than 2.5 times larger than in case of a non-rotating BH. The 
difference is smaller for larger torus masses. Considering 
that low-mass accretion tori rotating around a BH with 
intermediate values of the dimensionless angular momen- 
tum parameter (a ~ 0.6 — 0.8) may be typical products 
of mergers of comp act objects (e.g., iRuffert et al.l Il997t 
ISetiawan et al.|[2004[ ). we have demonstrated that all other 
GR effects are of moderate importance when calculating 
the amount of energy released in such systems by vv- 
annihilation. However, we point out that the spatial dis- 



tribution of the released energy exhibits very important 
differences between Newtonian and GR models: Close to 
the rotation axis and close to the bounding of the 'up' re- 
gion, GR effects can enhance the local energy deposition 
rate by a factor of ~ 10. Thus, a relativistic jet driven by 
the energy release in a region close to the stagnation sur- 
face (which will form in the vicinity of the boundary of 
the 'up' region) should receive a much larger energy input 
due to GR effects. The question whether this difference in 
the spatial distribution of the energy deposition is more 
favorable for producing ultrarelativistic jets can not be ad- 
dressed by the present work, but requires time-dependent 
hydrodynamic jet simulations including a detailed descrip- 
tion of i^i/- annihilation (see below). 

In our set of idealized models, a thin disk (sphere) is 
the most (least) favorable geometry to convert neutrino 
energy into e+e~-pairs by i/z/-annihilation. Our idealized 
models yield efficiencies of only a few tenths of a 

per cent. Larger efficiencies, of the order of several per- 
cent, have been obtained for our more realistic equilibrium 
models. The se findings a r e in t he ball-park of the results 
obtained by iJaroszvnskil (Il996l) ( f or mo dels with specific 
entropies > 8 and ISetiawan et all (|2004f ): for models with 
the largest a- viscosity) . 

IMiller et all (j2003D have pointed out that the mam con- 
tribution to the energy-momentum deposition rate comes 
from large polar angles (particularly from regions above the 
neutrinospheres) and not from regions near the symmetry 
axis. This behavior, however, results from their restriction 
to infinitely thin disk models, which are mathematically 
idealized cases. According to our studies, only such ideal- 
ized models show this strong enhancement of the annihi- 
lation rate towards the equatorial plane, while this effect 
is absent in the mor e realistic equilibri um torus models. 
Hence, this finding of IMiller et al.l (j2003l ) appears to be of 
little relevance for the amount of energy that is available 
to produce a GRB. In any realistic situation the accretion 
disk/torus will be inflated vertically, i.e. the baryon density 
will be highest in the equatorial plane and decrease in per- 
pendicular direction, the density scale height of the baryon 
distribution being larger for large disk masses and high disk 
temperatures. 

As demonstrated in lAlov et al.l (j2005[ ). who considered 
similar equilibrium tori as initial models for their hydro- 
dynamic simulations of relativistic jets from BH-torus sys- 
tems, time-dependent numerical simulations are needed in 
order to determine the feedback of z/j/-energy deposition on 
the torus structure, which, in its turn, can yield a num- 
ber of highly non-lin ear relativistic hydrod ynamic effects 
on the outflowing jet (|Alov fc Rezzollall2006( ). This also de- 
cides about which fraction of the energy deposited by vv- 
annihilation is finally useful for accelerating an ultrarela- 
tivistic fireball powering a GRB event. This fraction can 
be different from the annihilation efficiency defined 
in the present work. Our equilibrium models produce vv- 
annihilation rates E^y^'°" > lO^^ergs^^ for torus masses 
™tor > O.IMq. These figures agree within a factor of 2-4 
with the results of dynamical simulations (|Setiawan et al.l 
I2OO5I ) and they are abo ut one order of ma gnitude higher 
than the ones quoted in IJaroszvnskil (jl996D for his models 
with specific entropies per baryon in th e range of 8 to 10. 
The difference compared to I Jaroszvnskif s estimates is prob- 
ably due to the differences in the temperature distribution 



R. Birkl etal.: Neutrino pair anniiiilation near accreting black holes 



17 



on the neutrinosphere caused by our somewhat different 
way of constructing the accretion tori. The neutrino lumi- 
nosity is very sensitive to the neutrinospheric temperature, 
and a ^ 30% larger value of this temperature can account 
for a factor of more than 10 larger total i^z/-annihilation 
rate. Considering a certain accretion disk mass and its cor- 
responding neutrino luminosity and annihilation cfficicney 
means an upper bound of the physically likely situation, in 
which the neutrino luminosity decreases with time as the 
mass of the disk decreases. However, we remark that the 
time evolution of the ^'i/-annihilation rate and the annihila- 
tion efficiency can be non-monotonic and do es not necessar- 
ily de crease with time as demonstrated by ISetiawan et al.l 

When attempting to link our present results to GRB ob- 
servations, we must therefore be aware of the restrictions 
stated in the previous paragraph and in particular of the 
fact that any reliable estimate of the fraction of the lyD- 
annihilation energy that drives an ultrarelativistic fireball 
and the determination of the coUimation of such outfiow 
requi res hydrodynamic simulations (see also iJanka et al.l 
\20m . So we conclude with caution that the energy release 
by i/P-annihilation in some of our models (especially those 
having the most massive tori or large values of a) could be 
sufficient even to fuel the most distant and most powerful 
short GRB discovered so far ( GRB060121, « 3-105^erg; 
Ide Ugarte Postigo et al.|[2006l ) without the need of invoking 
either ultra-intense magnetic fields {B > 10^^ G) or a dif- 
ferent progenitor class. 

Finally, we point out that lAlov et al.l employed a simple 
fit with time-independent geometry to describe the energy 
deposition by j/j>-annihilation above the poles of a stellar- 
mass BH in their relativistic hydrodynamic jet simulations. 
Preferably, future numerical simulations should include a 
time-dependent treatment of the energy deposition by I'D- 
annihilation using the refined methods employed in this 
work. 
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Appendix A: Some technical details of the 
calculation of the annihilation rate 

Here we describe in detail the evaluation of Eq. ^ , which 
is a natural generaliza tion of the fo r mula for the energy 
component Q* given in iRuffert et all (I1997D . 

We express the neutrino 3-momentum p (analogous ex- 
pressions hold for the antineutrino quantities which are rep- 
resented by primed variables) in terms of the energy E = \p\ 
and the unit vector in the momentum direction n = p/E. 
The distribution function is defined as 



f = fit,x,E,n) 



dN 
JdV^dVp' 



where h is Planck's constant, and g is the statistical weight 
{g = 1 for neutrinos and antineutrinos). Using d^p — 
E^dEd^l and defining 9 as the angle between the directions 
of propagation of the neutrino and antineutrino, Eq. ^ in 



full detail reads 



4 TO^/l^ 

(Cl+C2),._ 



poo nOO 

/ dE dE' {p" +p'")E^E'^ 
Jo Jo 

(f dnf dn' (1 - cos 61)^ ^■ 

poo pOO 

C^.yru.ml dE dE'{p''+p"')E^E'^ 
Jo Jo 

f dnf dfl' {1~ cos 6) f^J^, 



where = {E,p), ctq = 1.76 • 10 '*'*cm^ is the weak inter- 
action cross section, mc the electron mass, and finally 

(Ci + C2),^,^ « 2.34, (Ci + C2),^,^ « 0.50, 

C3,.,p, » 1.06, C3,.,p, « -0.16, 

with X € {/i, r}. 

In order to perform the ray-tracing of the neutrino tra- 
jectories a definition of the base vectors spanning the local 
o bserver frame is req uired. Our choice is the same as that 
of iMiUer et all ((2003f) . i.e. 

(et") = ^ (1, 0, 0, 0) , (e,") = -jL= (0, 1, 0, 0) , 



9tt 
1 



(0,0,1,0), 



9tt 



9t<j3 ~ 9<P<t>9tt \ 9tt 



V-9r 



- — ,0,0,1 I , 



where the m etric coefficients are given in Eg. ([2l). No t e that 
iMiller et al.l used the same signature as lMisner et al.l (|1973f ) 
and that the selected base is orthonormal, i.e. ec, -e^ = rjajs, 
where r]ai3 is the Minkowski metric. 

Appendix B: Convergence Tests 

We performed a series of convergence tests to check the 
dependence of the energy component Q* of the annihila- 
tion rate 4- vector on the number N^ays of ray-tracing paths 
for neutrinos and antineutrinos. The results of this test 
are shown in Fig. IB. II The dots give the values of Q* at 
selected points for a representative BH-disk configuration. 
Successively increasing Arrays (gray filled circles connected 
by the black, solid line) we found that the value of Q* is 
converged to an accuracy better than ~ 3% for A'^rays > lO"*. 
When calculating Q* repeatedly for four selected values of 
iVrays the rcsults scatter statistically due to the random 
procedure used for picking the initial direction of the ray- 
tracing paths for the (anti)neutrinos. The corresponding 
relative error is shown in Fig. IB. II too. Since all the sim- 
ulations presented in this publication have been performed 



with Ni-s 



20000, 



the relative error of the calculated 
annihilation rates Q* is about 2%. This accuracy slightly 
varies depending on where the annihilation rate is calcu- 
lated. Analogous tests were also performed for the spatial 
components of the annihilation rate 4- vector, which show a 
similar convergence behaviour. 

A second series of convergence tests was computed, in 
order to check the dependence of the total annihilation rate 
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rays 



Fig. B.l. Convergence behaviour of the energy component 
Q* of the annihilation rate 4- vector for a M = IM0, 
a = black hole. The neutrinosphcre (for both and 
Da) is considered to be a non-rotating isothermal thin disk 
of temperature T = 10^^ K extending from r\a = 3.5-M to 
^out = 5.5A'/. The rate is computed at the point (r = '6 AM , 
= j) using different numbers -/Vj-ays of ray-tracing paths 
for neutrinos and antineutrinos. The horizontal line marks 



the asymptotic value of Q* as Nra 



on the grid resolution. In this context model REF, which 
was simulated on a computational grid of Nr = 110 and 
Ng = 100 points, was recalculated on a grid with twice as 
many grid points in each coordinate direction (i.e., Nr = 
220, and Ng = 200). Doubling the resolution the value of 
eI°p'°° defined in Eq. ^ changes from 6.19 • lO^^ergs"! 
(coarse grid) to 6.21 • 10'*^ergs~^ (finer grid), which is a 
small difference of 0.3%. 

Finally, we investigated the dependence of the results 
on the size of the grid, i.e. on the value of the outer grid 
radius ''g.out- Again using model REF, we increased the 
outer radius from its standard value out = 48.75 km to 



r'g^out = 270 km. The resulting difference in Ejl^ '°° is 2.5 • 
10^^ ergs~^, corresponding to a relative error of ~ 4%. This 
error is also representative for the other models (note that 
there is negligible energy deposition for r > rg ^^^.). 
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